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EXPERIMENTAL  FLUTTER  RESULTS  FOR  CORRUGATION-STIFFENED 

AND  UNSTIFFENED  PANELS 

By  Deene  J.  Weidman 
Langley  Research  Center 

SUMMARY 

Frhe  results  of  a  flutter  investigation  of  19  corrugation- stiffened  and  6  flat  unstiff¬ 
ened  rectangular  panels  are  presented^  These  results  show  that  boundary  conditions  in 
the  cross- stream  direction  have  an  irhportant  effect  on  panel  flutter.  It  is  also  shown 
that  the  flutter  of  square  orthotropic  panels  is  not  easily  predicted,  and  that,  at  present, 
a  flutter  envelope  appears  to  be  the  best  available  means  for  the  prediction  of  the  flutter 
boundaries  for  such  panels.^  In  addition,ia  method  for  analyzing  the  stiffness  of 
corrugation- stiffened  skin  panels  is  presented.J  A  "preflutter"  solution  is  shown  to 
approximate  accurately  the  exact  solution  for  simply  supported  panels  with  in-plane  ten¬ 
sion  loads,  but  both  solutions  differ  considerably  from  the  experimental  results. 

INTRODUCTION 

With  the  increased  speed  and  performance  of  high-speed  aircraft  and  missiles,  the 
problem  of  flutter  of  the  individual  skin  panels  of  such  vehicles  has  become  more  pro¬ 
nounced  and,  indeed,  has  even  governed  design  in  some  cases.  Many  skin  panels  of  these 
vehicles  are  constructed  of  an  exterior  sheet  stiffened  by  means  of  a  single  corrugated 
stiffening  sheet  attached  directly  to  the  skin.  These  panels  are  termed  "corrugation- 
stiffened"  panels  herein.  In  the  theoretical  investigation  of  panel  flutter,  several  basic 
approaches  have  long  been  available  (see  refs,  1  to  3),  and  recent  extensions  of  these 
analyses  (refs.  4  to  7)  indicate  a  continuing  interest  in  this  problem.  Experimental 
results,  however,  differ  considerably  from  most  analytical  results  particularly  for 
corrugation- stiffened  panels  (see,  for  example,  ref.  8),  and  the  design  of  panels  has,  to  a 
large  extent,  been  accomplished  with  the  use  of  an  experimental  flutter  envelope  such  as 
is  presented  in  reference  9. 

In  view  of  the  difficulty  in  theoretically  determining  the  flutter  characteristics  of 
corrugation- stiffened  panels,  it  appears  that  experimental  flutter  data  should  be  used  to 
determine  a  flutter  envelope  for  panel  design  until  a  proper  analytical  approach  can  be 
determined.  Isotropic  plate  theory  has  been  extended  to  orthotropic  plates  (if  the  ortho¬ 
tropic  stiffnesses  are  known)  as  in  references  8  and  9;  but  the  validity  of  such  an 


extension  is  questioned  in  reference  10.  However,  this  extension  defines  parameters 
that  allow  the  experimental  results  to  be  compared  in  a  logical  manner;  these  param¬ 
eters  are  used  herein.  The  orthotropic  stiffnesses  of  these  corrugation- stiffened  panels 
are  calculated  by  the  method  shown  in  appendix  A. 

In  this  report,  the  results  of  a  flutter  investigation  of  both  isotropic  and 
corrugation- stiffened  panels  are  presented.  In  the  tests  of  corrugation- stiffened  panels, 
both  aerodynamic  and  panel  stiffness  factors  were  varied,  as  well  as  the  boundary  condi¬ 
tions  on  the  ends  of  the  corrugations.  In  all  tests  reported  herein,  the  axis  of  the  corru¬ 
gations  was  always  perpendicular  to  the  direction  of  airflow,  so  that  the  boundary  condi¬ 
tions  on  the  ends  of  the  corrugations  are  termed  "cross-flow”  boundary  conditions.  The 
pressure  differential  (difference  between  the  cavity  pressure  and  the  free- stream  pres¬ 
sure)  was  also  varied  for  these  panels.  In  the  tests  of  isotropic  panels,  the  aerodynamic 
parameters  and  pressure  differential  were  varied,  as  well  as  the  length-width  ratio  and 
thickness  of  the  panel. 

The  experimental  results  are  compared  with  the  empirical  flutter  envelope  from 
reference  9  and  with  the  "exact"  preflutter  solution  in  which  the  exact  eigenvalue  equa¬ 
tion  for  the  simply  supported  case  is  solved  for  a  conservative  critical  flutter  value. 

This  preflutter  solution  is  compared  with  the  exact  solution  and  shown  to  be  very  accu¬ 
rate.  A  typical  flutter  boundary  for  free-edge  panels  from  reference  5  is  also  presented. 
The  experimental  data,  however,  differ  considerably  from  these  solutions,  and  these  dif- 
ferences  are  discussed. 

SYMBOLS 


The  units  used  for  the  physical  quantities  in  this  paper  are  given  both  in  the  U.S. 
Customary  Units  and  in  the  International  System  of  Units  (SI).  Factors  relating  the  two 
systems  are  given  in  reference  11,  and  those  used  in  the  present  investigation  are  pre¬ 
sented  in  appendix  B. 


Nxa2  ^_„_2  a2  ^xy 


2n2^2^_±£ 

b2  Dx 


dimensionless  enclosed  area  within  a  cor¬ 
rugation  cell  (see  eq.  (A8)) 

length  and  width,  respectively,  of  each  indi¬ 
vidual  flutter  panel  (see  sketch) 


2 


E  Young's  modulus  of  corrugated  panel  material 

g  =  d  -  r(l  -  cos  0o) 


r  +  r 

hi,h2  quantities  in  stiffness  expression  (A6)  (defined  after  eq.  (A6)) 

K  =  w  -  S  -  S 
2(r  +  r) 

Z  =  d  -  (r  +  r)(l  -  cos  Bq) 

M  Mach  number 


Mx  bending  moment  applied  normal  to  x-direction 

Mj,M2  bending  moments  applied  to  skin  material  and  corrugation  material, 
respectively 
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m  integer  appearing  in  preflutter  solution 

Nx,Ny  applied  in-plane  forces  for  orthotropic  plate  (positive  in  compression)  in 
X-  and  y- directions,  respectively 

n  integer  appearing  in  exact  flutter  solutions 

P  applied  load  in  stiffness  analysis 

q  dynamic  pressure 

r,r  radii  of  fillet  curves  on  corrugation  (see  fig.  11) 

S,S  dimensions  of  corrugated  panel  (see  fig.  11) 

Si,S2,S3  widths  of  equivalent  corrugation  (see  fig.  13) 

ti,t2  thicknesses  of  skin  material  and  corrugation  material,  respectively 

ui,U2,U3  deflections  of  equivalent  corrugation  (see  fig.  13) 

Wn  deflection  function 


pitch  of  corrugation 


Cartesian  coordinates  of  orthotropic  plate 


parameters  appearing  in  exact  flutter  solution 


/3  =  -  1 


length  of  portion  of  corrugation  (see  fig.  11) 
pressure  differential  across  panel  (positive  outward) 


slope  of  corrugation  (see  fig.  11) 


dynamic-pressure  parameter. 


experimental  dynamic-pressure  parameter 
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preflutter  dynamic-pressure  parameter 


|Li  mass  of  orthotropic  plate  per  unit  area 

u  Poisson's  ratio  (assumed  to  be  1/3  for  numerical  calculations) 

CO  circular  frequency 

PANELS,  APPARATUS,  AND  TEST  PROCEDURE 

Panels 

Flutter  results  for  two  types  of  flat  rectangular  panels  are  presented:  corrugation- 
stiffened  panels  and  unstiffened  panels.  These  two  types  are  discussed  separately. 

Corrugation- stiffened  panels.-  Two  series  of  corrugation- stiffened  panels  were 
tested.  The  pertinent  dimensions,  materials,  and  stiffnesses  (calculated  by  the  method 
of  appendix  A)  are  given  in  table  I.  All  these  panels  were  tested  with  the  airflow  direc¬ 
tion  normal  to  the  axis  of  the  cormgations.  The  panels  were  assumed  essentially  to  be 
simply  supported  in  the  flow  direction  (that  is,  at  the  leading  and  trailing  edges). 

The  nine  panels  of  the  first  series  had  the  ends  of  each  corrugation  essentially  free 
and  are  called  free- edge  panels  herein.  Details  of  these  panels  are  shown  in  figure  1. 


Spacer 


Figure  L-  Corrugation-stiffened  panels  with  ends  of  corrugations  free  along  longitudinal  edges.  All  dimensions  are  given  first  in  inches  and 

parenthetically  in  centimeters. 
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The  panels  were  24,0  inches  (61,0  cm)  square,  and  had  the  corrugations  ending  approxi¬ 
mately  0.25  inch  (0,63  cm)  from  the  supporting  edge,  with  only  the  thin  skin  carried  over 
the  support  (see  section  A- A).  These  panels  are  designated  by  the  letter  "U,"  "V,"  or 
"H"  (hat),  which  describes  the  cross-sectional  shape  of  a  single  corrugation. 

Details  of  the  ten  panels  of  the  second  series,  on  which  the  corrugations  were 
directly  attached  to  the  support  fixture  at  their  ends,  are  shown  in  figure  2.  The  attach¬ 
ment  was  accomplished  by  crushing  the  ends  of  the  corrugations  and  allowing  this  thick 
crushed  region  to  extend  over  the  panel  support.  These  panels  are  called  crushed-edge 
panels  herein,  and  were  37.3  inches  (94.7  cm)  long  and  24.2  inches  (61.5  cm)  wide  with 
S-shaped  corrugations.  (See  fig.  2.)  These  crushed-edge  panels  are  designated  by  the 
letter  "S." 

The  panels  designated  S  and  V  had  a  single  weld  line  between  the  corrugations, 
whereas  the  panels  designated  U  and  H  had  double  weld  lines. 


Figure  2.-  Corrugation-stiffened  panels  with  corrugations  crushed  aiong  iongitudinai  edges.  All  dimensions  are  given  first  in  inches  and 

parenthetical iy  in  centimeters. 
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Unstiffened  panels.-  Two  series  of  flat  unstiffened  rectangular  panels  were  tested. 
The  details  of  these  panels  are  shown  in  figure  3,  and  the  materials  and  dimensions  are 
listed  in  table  n.  The  four  panels  of  the  first  series  were  divided  into  several  bays  in 
the  cross-flow,  direction  and  were  termed  "multi-bay  panels."  The  two-edge  bays  were 
equal  in  width  but  narrower  than  the  remaining  bays  in  order  that  the  total  width  of  each 
panel  assembly  remain  24.2  inches  (61.5  cm).  The  length-width  ratios  of  these  multi-bay 
panels  were  approximately  6  and  10.  The  two  panels  of  the  second  series  were 
24.2  inches  (61,5  cm)  wide  and  were  termed  "single-bay  panels"  (with  length-width 
ratios  a/b  of  1  and  1.54). 


24.2 (61 ,5) 


a 


(a)  Multi-bay  panels. 


(b)  Single-bay  panels. 

Figure  3.-  Dimensions  of  flat  unstiffened  panels.  All  dimensions  are  given  first  in  inches  and  parenthetically  in  centimeters. 


Test  Setup 

The  panel  flutter  tests  discusse(i  herein  were  conducted  at  the  Langley  Unitary  Plan 
wind  tunnel  in  the  fixture  shown  in  figures  4  and  5  and  discussed  in  reference  12.  This 
test  fixture  allows  any  panel  up  to  approximately  3  feet  (0.91  m)  long  by  2  feet  (0.61  m) 
wide  to  be  mounted,  the  reversible  vertical  splitter  plate  holding  the  panel  33  inches 
(0.84  m)  from  the  glass  observation  wall.  (See  fig.  5.)  The  pressure  distributions  for 
this  test  fixture  are  shown  in  reference  13  and  are  noted  there  as  being  nonuniform,  par¬ 
ticularly  at  the  low  Mach  numbers.  Tests  were  run  in  the  low  Mach  number  test  section 
of  the  Langley  Unitary  Plan  wind  tunnel  where  dynamic  pressure  can  be  varied 
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continuously  from  130  to  2600  psf 
(6  to  124  kN/m2)  and  Mach  numbers 
from  1.57  to  2.87  are  easily  obtainable. 
The  stagnation  temperature  of  this  tun¬ 
nel  can  be  varied  slightly,  but  was  held 
at  125°  F  (325°  K)  for  most  of  the  tests 
reported  herein. 

The  pressure  loading  across  the 
panels  was  also  varied  in  these  tests. 

The  pressure  differential  is  considered 
to  be  the  difference  in  pressure  between 
the  cavity  behind  the  panel  and  a  static - 
pressure  orifice  located  9  inches 
(23  cm)  ahead  of  the  leading  edge  of  the 
panel.  A  positive  pressure  differential 
indicates  a  higher  pressure  in  the  panel 
cavity.  This  pressure  differential  can  be  regulated  or  held  constant  (as  the  test  requires) 
by  means  of  a  pump  and  servo  arrangement;  however,  it  is  not  completely  uniform  over 
the  entire  panel  surface.  (See  ref.  13.)  In  all  tests,  six  variable- reluctance-type 
deflectometers  and  an  oscilloscope  (with  associated  equipment)  were  used  to  monitor  the 
motion  of  the  panels.  Both  high-speed  (16-mm)  and  low-speed  (8-mm)  motion  pictures 
were  taken  to  observe  the  motion  of  the  panels  at  flutter. 

Test  Procedure 

Since  it  is  easier  to  vary  the  dynamic  pressure  than  to  change  the  Mach  number  in 
the  tunnel,  the  following  test  procedure  was  adopted:  For  any  given  panel,  a  low  Mach 
number  was  selected  and  either  a  low  dynamic  pressure  or  large  differential  pressure 
was  established.  Then,  holding  the  dynamic  pressure  and  Mach  number  constant,  the 
pressure  differential  across  the  panel  was  varied  in  an  attempt  to  initiate  flutter.  Flutter 
was  defined  as  the  combination  of  aerodynamic  forces  and  pressure  differential  that 
caused  the  panel  motion  to  become  regular  and  its  amplitude  to  exceed  a  predetermined 
arbitrary  magnitude  which  was  selected  to  exclude  forced  dynamic  motion  but  not  flutter 
(whether  or  not  destruction  of  the  panel  occurred).  If  flutter  did  not  occur,  the  dynamic 
pressure  was  increased  by  a  preselected  amount  and  the  pressure  differential  was  varied 
again.  This  procedure  was  continued  until  the  flutter  conditions  were  determined  for  the 
panel  at  that  Mach  number.  Then  the  Mach  number  was  increased  a  selected  amount  and 
the  procedure  was  repeated  again. 


L-59-3993.1 

Figure  4.-  Test  fixture  with  panel  mounted  in  test  section. 
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RESULTS  AND  DISCUSSION 


Corrugation-Stiffened  Panels 

Typical  results  of  tests  of  corrugation- stiffened  panels  are  shown  in  figure  6,  where 
the  variation  of  dynamic  pressure  and  pressure  differential  are  shown  for  several  cases. 
For  a  given  panel  and  Mach  number,  the  tests  proceeded  in  the  manner  indicated  by  the 
arrows.  Straight  lines  indicate  no  flutter  conditions  and  serrate  lines  indicate  flutter. 

The  lowest  dynamic  pressure  at  which  a  given  panel  fluttered  was  defined  as  the  critical 
dynamic  pressure  for  that  panel  at  that  Mach  number.  The  pressure  differential  across 
the  panel  at  the  initiation  of  flutter  was  used  to  define  the  actual  critical  flutter  point  as 
shown  in  figure  6.  However,  the  pressure  differential  is  only  approximate  since  the  flow 
conditions  over  the  panel  are  not  completely  uniform.  (See  ref.  13.)  Generally,  from 
plots  such  as  figure  6,  it  can  be  seen  that  an  increase  in  pressure  differential  (either 
positive  or  negative)  stopped  the  flutter  motion  of  the  panel. 

In  some  cases,  approximately  the  same  magnitude  of  either  positive  or  negative 
pressure  differential  was  needed  to  stop  flutter  (see  figs.  6(a)  and  6(b)).  In  other  cases, 
however,  a  somewhat  erratic  behavior  occurred,  as  shown  in  figure  6(c),  and,  for  these 
cases,  large  positive  pressure  differentials  were  generally  required  to  initiate  flutter. 
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(a)  Panel  V-2;  M  =  1.96.  (W  Panel  V-5;  M  =  1.84. 

Figure  6.-  Typical  test  data  for  corrugation-stiffened  panels. 
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The  reason  for  these  results 
could  be  initial  imperfections 
in  the  panels.  However,  this 
information  cannot  be  deter¬ 
mined  from  the  data  since*  no 
means  for  measuring  such 
effects  were  provided. 


The  results  of  the  flut¬ 
ter  tests  of  corrugation- 
stiffened  panels  are  presented 
in  table  III,  where  for  each 
panel  at  specific  Mach  num¬ 
bers  the  values  of  q  and 
Ap  for  the  ^'critical  flutter 
point”  are  presented.  These 
results  are  also  shown  plotted 
in  figure  7  in  terms  of  the 
critical  dynamic-pressure 
flutter  parameter 
2qa3 


where  X  = 


/3Dx 


and  a 


parameter  A/4^.  In  the  gen¬ 
eral  case  (for  panels  simply 
supported  on  all  edges),  A 
can  be  written  as 


-2  0  2  4  6  8 

I - r - \ - 1 - 1 - 1 


(c)  Panel  V-6;  M  =  1.84. 
Figure  6.-  Concluded. 


A  = 


Nxa2 

Dx 


2n^Tr^ 


b2  Dx 


(1) 


which  combines  a  length-width  ratio  term  and  an  in-plane  loading  term.  Herein,  a  and 
b  refer  to  the  length  and  width  of  the  panel,  respectively.  is  the  in-plane  com¬ 

pressive  load  in  the  airflow  direction,  is  the  bending  stiffness  in  the  airflow  direc¬ 
tion,  and  Dxy  is  the  twisting  stiffness  of  the  panel.  Since  the  in-plane  forces  could  not 
be  either  directly  varied  or  measured  in  this  investigation,  the  in-plane  forces  were 
assumed  to  be  zero  and  the  first  term  in  the  expression  for  A  (eq.  (1))  was  neglected 
when  the  test  results  were  plotted.  The  quantity  n  is  an  integer  associated  with  the 
number  of  half-waves  in  the  cross-flow  direction,  and  for  all  panels  reported  herein  n 
was  taken  to  be  1.  The  panel  stiffnesses  were  computed  by  the  method  shown  in  appen¬ 
dix  A  and  are  presented  in  table  I.  The  vertical  I-shaped  bars  in  figure  7  represent  the 
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scatter  of  data  measured  for  any  given  panel,  and  the  designation  by  each  bar  or  point 
refers  to  the  corrugation  shape  and  panel  number  given  in  table  I. 


A 

Tr2 


Figure  7.-  Results  of  flutter  tests  on  corrugation-stiffened  panels. 

The  large  scatter  of  flutter  results  for  a  single  panel  indicates  that  effects  other 
than  those  incorporated  in  X  and  A  which  are  important  in  panel  flutter  were  probably 
present  during  testing.  One  effect  is  the  pressure  differential  mentioned  previously. 

The  influence  of  boundary  conditions  at  the  ends  of  the  corrugations  (cross- stream  bound¬ 
ary  conditions)  are  also  seen  to  be  important,  since  the  crushed- edge  panels  (designated 
S)  fluttered  at  significantly  higher  values  of  the  dynamic-pressure  parameter  than  the 
free-edge  panels  (designated  U,  V,  and  H).  This  fact  is  further  verified  by  the  test 
results  from  another  investigation  of  corrugation-  stiffened  panels.  (See  ref.  8.)  In  ref¬ 
erence  8,  two  of  the  free-edge  panels  were  attached  to  the  support  fixture  at  the  corruga¬ 
tion  ends  by  angle  clips.  (See  fig.  8.)  The  two  panels  selected  for  this  adjustment  were 
U-2  and  H-1;  some  results  of  flutter  tests  on  these  panels  are  also  shown  in  reference  8. 
The  panel  U-2  (fixed)  should  flutter  (for  zero  thermal  stress)  at  a  dynamic  pressure  near 
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3400  psf  (163  kN/m2)  and  panel  H-1  (fixed)  should  flutter  at  a  dynamic  pressure  near 
2900  psf  (139  kN/m2)  both  at  a  Mach  number  of  3,0.  These  two  data  points  are  shown  in 
figure  7  as  circled  points.  As  can  be  seen  from  these  two  data  points,  attaching  the  ends 
of  the  corrugations  to  the  support  fixture  caused  significantly  higher  dynamic  pressures 
for  initiation  of  flutter  than  for  identical  free- edge  panels.  Therefore,  the  attachment  of 
the  ends  of  the  corrugation  greatly  increases  the  flutter  resistance  of  a  corrugation- 
stiffened  panel.  (This  result  is  also  noted  in  ref.  8.) 


Section  B— B 


Spacer- 


Ang I e  clip: 
-.75X.75X .030 
(1.91x1 .91 X. 076) 


r777777////////i 


Section  A— A 


Figure  8.-  Corrugation-stiffened  free-edge  panels  modified  to  approximate  clamped  edges  (from  ref.  81. 


This  effect  can  be  explained  by  assuming  that  the  panels  effectively  end  where  the 
corrugations  either  end  or  are  crushed.  The  remaining  portion  of  the  panel  must  be  con¬ 
sidered  as  a  flexible  support.  With  this  assumption,  the  free-edge  panels,  the  modified 
free-edge  panels,  and  the  crushed-edge  panels  could  be  idealized  as  shown  in  figure  9. 
The  three  boundary  conditions  investigated  can  be  considered  to  be  special  cases  of  an 
orthotropic  plate  supported  on  a  general  elastic- edge  support.  The  transverse  flexibility 
of  the  panel  support  determines  the  deflectional  spring  constant  for  the  general  case, 
whereas  the  rotational  flexibility  of  the  panel  support  determines  the  torsional  spring 
constant. 
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Figure  9.-  Boundary-condition  idealization. 


From  these  idealizations,  it  can  be 
seen  that  the  free- edge  panel  supports 
correspond  to  the  general  case  with  a 
relatively  small  amount  of  transverse 
stiffness  and  negligible  rotational  stiff¬ 
ness.  The  modification  of  the  free-edge 
panel  supports  (by  addition  of  an  angle 
clip)  supplies  a  rotational  stiffness  and 
increases  the  transverse  stiffness  some¬ 
what.  The  crushed- edge  panel  support 
supplies  not  only  rotational  stiffness  but 
also  an  appreciable  amount  of  transverse 
stiffness  as  well.  Each  stiffness  increase 
causes  an  increase  in  the  critical  dynamic 
pressure  for  flutter.  All  these  conclu¬ 
sions  are  borne  out  by  the  experimental 
results  shown  in  figure  7.  A  similar  con¬ 


clusion  was  noted  in  reference  14,  where  improved  agreement  between  theory  and  experi¬ 
ment  was  obtained  with  the  use  of  an  estimated  value  of  this  spring  constant.  The  effect 
of  these  edge  conditions  is  discussed  in  the  next  section. 


The  crushed-edge  panels  also  fluttered  at  higher  values  of  the  dynamic-pressure 
parameter  X  than  the  free-edge  panels.  The  mode  of  failure  for  these  two  groups  of 
panels  was  entirely  different.  The  crushed- edge  panels  failed  quickly  and  catastrophi¬ 
cally,  and  the  downstream  ends  of  the  panels  were  destroyed.  The  free-edge  panels,  on 
the  other  hand,  fluttered  for  some  time  at  several  amplitude  levels  before  small  cracks 
formed  in  the  skin  at  the  ends  of  the  corrugations.  These  panels  still  maintained  their 
integrity  even  after  cracks  appeared.  A  similar  dependency  of  mode  of  failure  on  the 
flexibility  of  edge  supports  were  noted  in  reference  8. 


Unstiffened  Panels 

The  results  of  the  flutter  tests  on  flat  unstiffened  single-bay  and  multi-bay  panels 
are  presented  in  table  IV.  The  panel  numbers  are  the  panel  designations  given  in 
table  n.  For  each  Mach  number  used,  the  dynamic  pressure  and  pressure  differential 
on  the  panel  at  the  flutter  are  listed  for  the  conditions  discussed  under  the  column  heading 
"Flutter  Motion  and  Comments." 

These  test  results  are  also  shown  plotted  in  figure  10  in  terms  of  X^/^  and  A/4^. 
In  plotting  the  data,  the  in-plane  forces  were  assumed  to  be  zero,  as  was  done  with  the 
corrugation- stiffened  panels.  The  vertical  I-shaped  bars  in  figure  10  again  represent 
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scatter  in  the  results,  and  indicates  that  other  effects  (pressure  differential,  for  example) 
seem  to  be  important  for  flat  unstiffened  panels. 


Figure  10.-  Results  of  flutter  tests  on  flat  unstiffened  panels. 


COMPARISON  OF  TEST  RESULTS  WITH  THEORY 

As  a  simple  means  for  comparing  experimental  results  with  theory,  a  preflutter 
solution  can  be  written  that  is  both  easy  to  compute  and  accurate  in  the  region  of 
corrugation- stiffened  panels.  The  "exact"  solution  from  reference  2  assumes  small 
deflections  and  thin  -plate  theory  with  static  aerodynamic  forces  for  a  simply  supported 
panel.  The  differential  equation  for  the  deflection  (eq.  (5)  of  ref.  2)  then  becomes 

+  AWn  +  XW^  -  BWn  =  0  (2) 
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When  this  equation  is  solved  exactly,  it  yields  the  eigenvalue  equation  (eq.  (9)  of  ref.  2) 

(e  2  sin  6  sinh  e  =  5(cosh  e  cos  6  -  cosh  2a)  (3) 


where 


and  Q!  is  a  parameter  that  is  varied  until  equation  (3)  is  satisfied  for  any  given  values 
of  A  and  X.  The  definition  of  the  quantities  in  equation  (2)  for  orthotropic  plates 
(from  ref.  9)  are  given  in  the  notation  of  the  present  paper  as 


A  = 


Nx3' 


.2„v4^ 

Dx 


B 


2^  2  2  9  3-^  4  4  a^ 

=  Ui:^CO^  +  n^TT^  — - n'^IT^  —  — 


Px 


Dx 


b^ 


_  2qa^ 

/3Px 


By  following  the  method  of  Movchan  (ref,  15),  a  solution  of  the  "exact"  eigenvalue  equa¬ 
tion  (eq.  (3))  which  yields  conservative  estimates  of  the  flutter  boundary  can  be  shown  to 
exist.  Equation  (3)  is  thus  seen  to  be  satisfied  if  6  =  ±2m'ir  and  e  =  2q!.  Substituting 
these  values  into  equations  (4)  yields  the  preflutter  solution  for  the  critical  flutter  param¬ 
eter  as 


Xp  =  8  (m2,r2  + 

For  panels  with  the  axis  of  the  corrugations  normal  to  the  direction  of  flow,  the  twisting 
stiffness  Dxy  is  much  larger  than  the  bending  stiffness  Dx,  and,  thus,  for  small 
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compressive  forces  Nx,  negative  values  of  A  result.  (See  eqs.  (5),)  In  this  case, 

Xp  has  a  minimum  value  when  m  =  1  (this  value  of  m  corresponds  to  the  coalescence 
of  the  first  two  mode  shapes  only  and  may  be  somewhat  in  error).  When  m  =  1,  equa¬ 
tion  (6)  becomes 


X 


P 


16.878 


(7) 


This  preflutter  solution  is  compared  in  the  following  table  with  the  exact  solution  of  the 
eigenvalue  equation  (eq.  (3)): 


=4  l> 

^ex 

-0 

337.56 

341 

-10 

1  263,0 

1  314 

-20 

2  480.6 

2  602 

-100 

18  932 

19  390 

-1  000 

540  120 

542  220 

-10  000 

16  898  000 

16  904  000 

From  this  comparison,  it  can  be  seen  that  the  preflutter  solution  does  agree  very  closely 
with  the  "exact"’  solution  for  large  negative  values  of  A  for  simply  supported  panels. 
This  preflutter  solution  may  now  be  compared  with  the  experimental  results  to  determine 
the  actual  applicability  of  such  an  approach. 

Corrugation- Stiffened  Panels 

The  experimental  flutter  results  for  the  corrugation- stiffened  panels  presented  in 
table  III  are  compared  with  the  exact  preflutter  solution  in  figure  7.  From  the  figure  it 
is  seen  that  the  agreement  between  the  experimental  results  and  the  exact  preflutter  solu¬ 
tion  for  simply  supported  panels  is  extremely  poor,  the  agreement  becoming  even  worse 
for  large  negative  values  of  A.  Since  it  appears  from  the  test  results  that  boundary  con¬ 
ditions  in  the  cross-flow  direction  can  be  important,  an  additional  theoretical  curve  (from 
ref.  5)  has  been  added  to  figure  7.  This  curve  is  the  result  of  an  approximate  two-mode 
Galerkin  solution  for  free- edge  panels  with  stiffnesses  roughly  equal  to  those  of  the 
corrugation- stiffened  panels  tested.  The  experimental  data  are  seen  to  be  between  the 
curve  for  the  preflutter  solution  for  the  simply  supported  panels  and  the  curve  for  the 
Galerkin  solution  of  the  free-edge  panels.  This  fact  suggests  that  the  test  panels  were 
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somewhere  between  completely  free  and  simply  supported  with  respect  to  deflection  of 
the  corrugation  ends. 

The  experimental  flutter  envelope  from  reference  9  is  shown  in  figure  7;  this  curve 
also  falls  between  the  two  theoretical  curves.  Even  though  use  of  the  flutter  envelope  of 
reference  9  as  a  design  criterion  could  result  in  unnecessary  weight  penalties  in  some 
cases  with  a  possibility  of  being  unconservative  on  occasion,  this  type  of  curve  appears 
to  be  the  best  available  means  for  the  prediction  of  the  flutter  boundary  for  corrugation- 
stiffened  panels  at  present.  It  may  be  noted  from  figure  7  that  the  crushed- edge  panels 
generally  flutter  as  a  group  at  values  of  above  those  of  the  flutter  envelope, 

whereas  all  the  free-edge  panels  fluttered  at  values  of  X^/^  below  the  flutter  envelope. 
This  result  indicates  that  the  experimental  envelope  is  a  reasonable  curve  for  elastically 
supported  panels.  It  should  also  be  pointed  out  that  some  of  the  difference  between  the 
curve  for  simply  supported  panels  and  the  experimental  results  could  also  be  due  to  the 
fact  that  the  actual  transverse  shear  stiffnesses  of  these  panels  are  not  infinite  as  is 
assumed  in  the  theoretical  calculations.  The  influence  of  these  shearing  stiffnesses 
could  be  appreciable  and  this  problem  is  discussed  in  reference  5. 

Unstiffened  Panels 

The  experimental  flutter  results  for  the  flat  unstiffened  panels  presented  in 
table  IV  are  compared  in  figure  10  with  the  preflutter  solution  presented  previously. 

This  preflutter  solution  is  only  in  fair  agreement  with  the  experimental  flutter  results 
and  is  unconservative  for  larger  negative  values  of  a/t^.  The  experimental  envelope 
from  reference  9  is  also  shown  plotted  in  figure  10,  and  agreement  between  this  theoret¬ 
ical  curve  and  the  experimental  data  is  fairly  good,  particularly  at  larger  values  of  a/it^ 
This  curve  is  only  slightly  unconservative  in  this  region  and  can  be  considered  adequate 
throughout  the  entire  region  shown. 


CONCLUSIONS  , 

Upon  investigation  of  the  experimental  results  obtained  from  tests  Of  19  corrugation 
stiffened  panels  and  6  flat  unstiffened  panels,  several  important  facts  appear  evident:  ^ 

1.  A  flutter  envelope  drawn  to  encompass  all  experimental  results  available  at  any 
given  time  seems  to  be,  at  present,  the  best  means  for  predicting  the  flutter  boundary  for 
corrugation- stiffened  panels.  Theoretical  methods  currently  available  are  not  adequate 
to  predict  the  flutter  characteristics  of  a  general  panel.  The  flutter  envelope  from  NASA 
TN  D-451  agrees  with  the  test  results  for  isotropic  panels  fairly  well.  The  agreement 
for  crushed-edge  orthotropic  panels  is  fair,  but  the  agreement  for  free-e(^e  orthotropic 
panels  is  poor.  The  reason  for  this  poor  agreement  is  that  panels  with  free  edges  were 
not  used  to  define  the  flutter  envelope. 
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2.  The  mode  of  failure  for  panels  with  free- edge  corrugations  is  vastly  different 
than  the  mode  of  failure  for  the  crushed-edge  panels.  The  free-edge  panels  started  to 
flutter  at  much  lower  dynamic  pressures  and  continued  to  flutter  for  a  considerable 
length  of  time  before  cracks  appeared  in  the  skin;  even  then  panel  integrity  was  main¬ 
tained.  The  crushed-edge  panels,  however,  withstood  a  much  higher  dynamic  pressure 
before  flutter,  but  when  flutter  occurred  it  was  quick  and  catastrophic. 

3.  The  flutter  results  are  not  functions  only  of  the  flutter  parameter  and  the 
in-plane  loading  factor  but  are  also  greatly  affected  by  several  other  effects  that  are  not 
incorporated  in  these  parameters.  Several  of  the  effects  that  should  be  included  are 

(a)  The  influence  of  pressure  differential  across  the  panel, 

(b)  The  influence  of  boundary  conditions  in  the  cross-flow  direction. 

The  experimental  results  show  that  changing  the  boundary  conditions  in  the 
cross-flow  direction  can  change  the  flutter  dynamic  pressure  by  a  large  fac¬ 
tor,  and  the  influence  of  these  boundary  conditions  must  be  investigated  v 

(c)  The  influence  of  using  a  finite  shear  stiffness  for  any  given  panel. 

The  corrugation- stiffened  panels  are  far  from  infinitely  stiff  when  resisting 
a  shearing  load,  and  thus  the  influence  of  this  shear  stiffness  might  be  of 
importance  for  this  type  of  flutter  panel. 

Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Langley  Station,  Hampton,  Va.,  October  6,  1965. 
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STIFFNESSES  OF  FLAT  PANELS  REINFORCED  BY  A  CORRUGATED  SHEET 


Since  it  is  a  common  practice  in  the  aircraft  industry  to  construct  skin  panels  of  a 
corrugation- sheet  combination,  the  flutter  characteristics  of  such  panels  must  be  inves¬ 
tigated.  These  panels  may  be  considered  as  anisotropic  plates  and  panel  flutter  theory 
may  be  applied;  however,  the  stiffnesses  of  the  panels  must  be  determined  first.  The 
method  of  approach  used  herein  is  essentially  the  same  as  the  method  presented  in  ref¬ 
erence  16.  Several  differences  are  noted  where  they  occur,  and  all  expressions  are 
written  in  an  easily  usable  form. 

A  schematic  drawing  of  a  single  cell  of  a  general  panel  composed  of  a  skin  sheet 
and  a  single  corrugated- stiff ener  sheet  is  shown  in  figure  11,  This  general  shape  can  be 
specialized  to  a  square  (hat)  corrugation,  a  V-corrugation,  or  a  U-corrugation  by  proper 
choice  of  the  parameters.  The  heavy  vertical  marks  between  skin  and  corrugation  indi¬ 
cate  either  weld  or  rivet  junctions  joining  the  two  parts. 


Bending  Stiffness  Dx 

One  approach  to  the  problem  of  computing  the  bending  stiffness  Dx  is  to  separate 
the  skin  sheet  and  the  corrugation  sheet  and  analyze  each  one  individually,  and  then  equate 
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their  slopes  and  deflections 
at  the  weld  line  so  that  they 
can  be  joined  again.  The 
skin  and  corrugation  are  sub¬ 
jected  to  the  loads  and 
moments  shown  in  figure  12. 

The  load  P  is  applied  so 
that  the  horizontal  displace¬ 
ment  of  the  end  of  the  corru¬ 
gation  (caused  by  the  applied 
moment  M2)  may  be  equated 
to  the  horizontal  displace¬ 
ment  of  the  cover  sheet.  By 
applying  Castigliano^s 
theorem  as  in  reference  17 
and  by  writing  the  strain 
energy  of  bending  for  both  the 
corrugation  sheet  and  the 
cover  sheet,  the  displacements  and 
slopes  of  both  sheets  can  be  deter¬ 
mined  at  the  weld  line.  Equating 
these  expressions  yields  two  equa¬ 
tions  involving  P,  Mj,  and  M2, 
and  these  equations  can  be  used  to 
solve  for  the  ratio  of  M^  to  M2. 
The  final  equation  needed  is  found 
by  assuming  that  the  stiffness  of  a 
panel  is  proportional  to  the 
moment  it  will  carry,  or 


Figure  12.-  Loads  and  moments  used  in  analysis  of  stiffnesses  of  a  general  corrugation. 


Dx  Ml  +  M2 


(Al) 


Figure  13.-  Idealization  of  corrugation  to  an  equivalent  corrugation. 


where  Dx/I^s  ^^e  ratio  of  the  stiffness  of  the  combination  of  the  skin  and  corrugation 
to  the  stiffness  of  the  skin  alone.  Thus  the  stiffness  is  determined  for  the  bent  portion  of 
the  corrugation  between  welds.  The  stiffness  of  the  short  section  of  "doubled"  plate 
between  the  welds  of  adjacent  corrugations  must  also  be  taken  into  account. 


There  are  several  ways  of  doing  so.  The  method  selected  herein  replaces  the  two 
sections  of  corrugated  panel  by  a  single  panel  section  with  an  equivalent  stiffness  so  that 
the  deflection  of  the  actual  panel  and  the  equivalent  panel  are  the  same  under  a  moment 
applied  at  the  tip,  as  shown  in  figure  13.  Thus, 
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/  2  2 
■^'^2  =  -^ls2  +2siS2  +  Sj  ^ 


=  Ui  +  u 


1  +  >^2 


Now  the  equivalent  stiffness  Dg  can  be  computed  if  the  component  stiffnesses  are 
known.  The  stiffness  of  the  bent  corrugation  Dg  has  been  calculated  from  equation  (A2) 
as  Dx-  The  stiffness  of  the  doubled  sheet  is  assumed  to  be  the  stiffness  of  a  plate 
of  thickness  tj  +  t2.  Thus,  from  equation  (A2),  the  equivalent  stiffness  ratio  Dx/Ds  in 
the  x-direction  is 


\h(  h  f 


w  /LDs\ti+t2/  J 

where  the  stiffness  of  the  corrugation  cell  per  unit  width  D^/d, 


w  -  S  t.  2ji 


2d  fe  -  c2 


where 
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in  00  d  ^Vd/  \d  dS  ° 


t;\3/30 


00  cos  00  sin-0o 


-  sin  20, 


with  L  =  1  -  +  lj(l  -  cos  0o] 


I  =  1  -  I(l  -  cos  0o). 
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It  may  be  noted  that  either  0q  of  S  need  be  given  but  not  both.  They  are  related  by 
the  following  expression: 

S  =  w  -  S  -  2  sin  eo(r  +  7}  -  2l  cot  9q 


or 


(A5) 


where  H  =  — ^  -  1  and  K  = 

r  +  r  2(r  +  r) 

For  panels  with  a  single  weld  line  between  corrugations  (y  =  s/2),  the  equivalent 
panel  stiffness  is  given  by  equation  (A4),  This  expression  is  equivalent  to  the  corre¬ 
sponding  expression  (B19)  from  reference  16.  For  panels  with  a  double  weld  line  between 
corrugations,  equation  (A3)  must  be  used,  and  this  expression  differs  slightly  from  the 
corresponding  expression  (B33)  of  reference  16,  where  an  equal  change  in  slope  between 
the  two  ends  of  the  equivalent  section  was  used  as  a  criterion  instead  of  equal  deflection 
of  the  tip  of  the  equivalent  section  as  was  used  herein. 

Bending  Stiffness  Dy 

The  bending  stiffness  of  these  panels  in  the  direction  of  the  corrugations  Dy  can 
be  determined  by  calculating  the  moment  of  inertia  per  unit  width  of  a  panel  with  the 
cross  section  given  in  figure  11.  However,  it  has  been  shown  (ref.  4)  that  the  stiffness  of 
a  panel  perpendicular  to  the  airflow  (characterized  by  the  parameter  B)  has  little  direct 
effect  on  the  problem  of  panel  flutter.  All  the  corrugation- stiffened  panels  tested  herein 
had  their  corrugations  perpendicular  to  the  airflow,  and,  therefore,  the  stiffnesses  of 
these  corrugated  panels  were  not  computed  in  the  direction  of  the  corrugations.  If,  how¬ 
ever,  the  corrugations  were  alined  with  the  airflow,  the  stiffness  Dy  would  be  important 
and  must  be  computed.  Therefore,  this  stiffness  expression  from  reference  16  is 
included  herein  for  completeness  and  is  as  follows; 
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Several  small  terms  involving  were  added  to  the  expression  from  reference  16  in 

order  that  equation  (A6)  could  be  consistent  and  include  all  terms  involving  !£  . 

Twisting  Stiffness  D^y 

The  twisting  stiffness  of  corrugation- stiffened  panels  can  be  computed  by  con¬ 
sidering  the  torsion-twist  relationship  for  each  corrugation  cell.  Thus,  in  reference  16, 
it  is  shown 


12Ao^(l  -  v) 

!i!W-  -  1  +  ^  +  2  ^  e 
^2pi\d  d+  d  ^  t2  , 


where  Aq  is  1/^2  times  the  averse  area  enclosed  by  the  inside  and  outside  of  a  cor 
rugation  cell.  This  nondimensional  area  can  be  written  as 
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Ao  = 


§  ,  Ir 

d  d 


+  2 


2Ad2 


+  cot  do  + 


2 


(A8) 


The  twisting  stiffness  expression  (eq.  (A7))  differs  from  the  corresponding  expression  in 
reference  9  (eq.  (11))  by  a  factor  of  2;  the  present  result  is  smaller.  This  difference, 
however,  is  insignificant  when  computing  flutter  values  for  corrugation- stiffened  panels 
since  the  flutter  boundary  curve  of  reference  9  is  very  flat  in  this  region.  This  fact  has 
been  noted  before.  (See  ref,  10.) 
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CONVERSION  OF  U.S.  CUSTOMARY  UNITS  TO  SI  UNITS 

The  International  System  of  Units  (SI)  was  adapted  by  the  Eleventh  General 
Conference  on  Weights  and  Measures,  Paris,  October  1960,  in  Resolution  No.  12  (ref.  11). 
Conversion  factors  for  the  units  used  herein  are  given  in  the  following  table: 


Physical  quantity 

U.S. 

Customary 

Units 

Conversion 

factor 

(*) 

SI  Units 

Length . 

Moment . 

Pressure  .... 

Temperature  .  . 

in. 

in-lbf 

lbf/ft2 

(OF  +  459.67) 

0.0254 

0.113 

47.88 

5/9 

meters  (m) 
meter-newton  (m-N) 

newtons/meter2  (N/m2) 

degrees  Kelvin  (OK) 

^Multiply  value  given  in  U.S.  Customary  Unit  by  conversion  factor  to  obtain 


equivalent  value  in  SI  unit. 

Prefixes  to  indicate  multiple  of  units  are  as  follows: 


Prefix 

Multiple 

kilo  (k) 
centi  (c) 

10^ 

10-2 
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TABLE  I.-  DIMENSIONS  AND  STIFFNESSES  OF  CORRUGATION- STIFFENED  PANELS 


TABLE  n.-  DIMENSIONS  OF  FLAT  UNSTIFFENED  FLUTTER  PANELS 


Panel 


Number  of 
bays 
(width  b) 


94.49 

6.00 

15.24 

94.49 

6.00 

15.24 

94.49 

3.62 

9.19 

94.49 

3.62 

9.19 

Multi-bay  panels 


160 


Single-bay  panels 

.26  0,1250  0,318 
.26  .0190  .048 


Material 


2024 -T3  aluminum  alloy 
2024- T3  aluminum  alloy 
2024-T3  aluminum  alloy 
2024-T3  aluminum  alloy 


2024 -T3  aluminum  alloy 
17-7  PH  stainless  steel 
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TABLE  in.-  FLUTTER  DATA  FOR  CORRUGATION-STIFFENED  PANELS 


From  reference  8. 


TABLE  IV.-  FLUTTER  DATA  FOR  FLAT  UNSTIFFENED  PANELS 


Panel 

Dynamic  pressure, 

q 

Pressure  differential, 

Ap 

Flutter  motion  and  comments 

psf 

kN/m2 

psf 

kN/m2 

Multi-bay  panels 

1 

2.10 

318 

1.523 

14 

0.067 

2.10 

366 

1.752 

12 

.057 

2.10 

417 

1.997 

10 

.048 

2.10 

467 

2.236 

6 

.029 

2.10 

517 

2.475 

4 

.019 

2.10 

566 

2.710 

1 

.005 

2.10 

244 

1.168 

9 

.043 

Increase  stagnation  temperature  50°  F  (28°  K) 

2.10 

147 

.704 

6 

.029 

Increase  stagnation  temperature  50°  F  (28°  K) 

1.86 

.780 

2 

.010 

1.86 

1.700 

4 

.019 

1.86 

2.236 

6 

.029 

1.86 

2.657 

12 

,057 

1.57 

.871 

9 

.043 

2 

2.10 

2.681 

17 

.081 

2.10 

2.801 

37 

.177 

No  flutter 

2.10 

2.394 

10 

.048 

2.10 

3.830 

10 

.048 

Large  amplitude  motion  initiated 

2.10 

6.004 

70 

.316 

2.10 

2.159 

«0 

2.10 

709 

3.395 

*0 

Large  amplitude  motion  initiated 

1.70 

822 

3.936 

«0 

1.70 

1283 

6.143 

~0 

Large  amplitude  motion  initiated 

3 

2.10 

1020 

4.884 

160 

.766 

2.10 

884 

4.233 

400 

1.915 

No  flutter 

4 

2.10 

1013 

4.850 

190 

.910 

Large  amplitude  motion  initiated 

Single-bay  panels 

5 

1.87 

1580 

7.565 

^0 

^0 

1.57 

1867 

8.939 

=^0 

^0 

6 

1.87 

1400 

6.703 

~0 

^0 

1.87 

1790 

8.571 

-0 

~0 

Failure  of  panel 
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NASA- Langley,  1966  1-1541 


”The  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomerta  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof.” 

— National  Aeronautics  and  Space  Act  of  1958 


NASA  SCIENTIFIC  AND  TECHNICAL  PUBLICATIONS 


TECHNICAL  REPORTS:  Scientific  and  technical  information  considered 
important,  complete,  and  a  lasting  contribution  to  existing  knowledge. 

TECHNICAL  NOTES:  Information  less  broad  in  scope  but  nevertheless 
of  importance  as  a  contribution  to  existing  knowledge. 

TECHNICAL  MEMORANDUMS:  Information  receiving  limited  distri¬ 
bution  because  of  preliminary  data,  security  classification,  or  other  reasons. 

CONTRACTOR  REPORTS:  Technical  information  generated  in  con¬ 
nection  with  a  NASA  contract  or  grant  and  released  under  NASA  auspices. 

TECHNICAL  TRANSLATIONS:  Information  published  in  a  foreign 
language  considered  to  merit  NASA  distribution  in  English. 

TECHNICAL  REPRINTS:  Information  derived  from  NASA  activities 
and  initially  published  in  the  form  of  journal  articles, 

SPECIAL  PUBLICATIONS:  Information  derived  from  or  of  value  to 
NASA  activities  but  not  necessarily  reporting  the  results  of  individual 
NASA-programmed  scientific  efforts.  Publications  include  conference 
proceedings,  monographs,  data  compilations,  handbooks,  sourcebooks, 
and  special  bibliographies. 


Details  on  the  availability  of  these  publications  may  be  obtained  from: 


SCIENTIFIC  AND  TECHNICAL  INFORMATION  DIVISION 

NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 

Washington,  D.C.  20546 


